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SUMMARY 


An  exploratory  study  is  made  in  connection  with  the  problem  of  pre¬ 
dicting  properties  of  the  turbulent  base  flow  region  of  a  body  at  hypersonic 
speeds.  An  attempt  is  made  to  calculate  the  gross  features  of  the  base  flow 
using  a  non-similar  mixing  model  including  the  Chapman-Korst  recompression 
condition.  The  expression  for  the  eddy  viscosity  is  based  on  a  model  which 
allows  the  turbulent  equations  to  be  transformed  into  incompressible  laminar 
form.  The  empirical  factor  in  the  eddy  viscosity  expression  is  evaluated 
from  data  obtained  from  the  near  field  (in  the  linear  growth  region)  of  jets 
exhausting  into  a  quiescent  region.  The  same  expression  for  the  eddy  vis¬ 
cosity  is  then  used  in  an  attempt  to  estimate  the  growth  of  non-similar 
turbulent  mixing  layer  using  both  turbulent  and  laminar  initial  profiles. 

The  correlation  of  turbulent  jet  mixing  data  for  speeds  up  to  Mach  3 
shows  that  the  eddy  viscosity  is  »  very  strong  function  of  Mach  number  (or 
of  density  ratio  across  the  jet,  since  the  data  are  for  adiabatic  flow). 

If  the  eddy  viscosity  dependence  on  density  ratio  persists  for  dens  ity 
ratios  typical  of  re-entry  conditions,  then  it  is  found  that  at  these 
conditions  the  growth  rate  of  the  turbulent  shear  layer  is  orders  of  magni¬ 
tude  slower  approaching  the  rate  of  growth  of  the  laminar  mixing  layer 
under  the  same  conditions.  This  is  in  maT ked  contrast  to  the  situation 
at  low  supersonic  speeds,  where  the  turbulent  mixing  layer  dividing  stream¬ 
line  velocity  approaches  the  similar  value  well  before  recompression  occurs. 
If  this  speculation  is  correct,  then  the  stagnation  enthalpy  at  recompress¬ 
ion  in  the  hypersonic  turbulent  base  flow  should  be  well  below  the  total 
enthalpy  value,  since  the  enthalpy  build-up  on  the  dividing  streamline 
will  be  similarly  suppressed. 


This  research  is  a  part  of  PROJECT  DEFENDER,  sponsored  by  the  Advanced 
Research  Projects  Agency,  Department  of  Defense,  under  ARPA  Order  No. 
234-62,  monitored  by  Air  Force  Ballistic  Systems  Division  under  Contract 
No.  AF  04(694) -570. 
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1.  INTRODUCTION 


Predictions  of  observables  in  the  wakes  of  slender  re-entry  vehicles 
have  been  found  to  depend  crucially  on  the  flow  properties  assumed  at 
the  neck  region.  For  example,  predicted  electron  densities  and  radar 
cross-sections  of  wakes  may  vary  by  orders  of  magnitude,  depending  on 
the  assumed  value  of  the  centerline  neck  enthalpy  ratio  h/H^.^  ^  Since 
properties  at  the  wake  neck  are  determined  by  the  flow  processes  in  the 
body  boundary  layer  and  base  flow  mixing  regions,  it  is  clear  that 
accurate  methods  for  flow  calculations  in  these  regions  must  be  developed. 

The  primary  objective  of  the  present  paper  is  to  devise  a  procedure 

for  calculating  the  gross  properties  of  the  hypersonic  turbulent  base 

flow  (e.g.  base  pressure,  neck  enthalpy,  etc.).  The  problem  is  approached 

from  the  point  of  view  of  the  generalized  Chapman-Korst  base  flow  model, 

(2  3) 

as  was  done  in  the  laminar  case.  ’  For  laminar  flow  this  model  is 

thought  to  be  reasonable  for  moderate  supersonic  Mach  numoers,  but  is 

probably  not  adequate  for  re-entry  conditions.  Recently,  Reeves  and 
(4)  (5) 

Lees,  Webb,  Go! ik  and  Lees,  and  others  have  attempted  to  devise 
base  flow  models  more  applicable  to  re-entry  conditions.  These  methods 
are  not  yet  well  established  and  still  may  not  contain  essential  features 
even  for  the  laminar  case.  For  turbulent  flow  the  difficulty  is  compounded 
by  uncertainties  related  to  the  "eddy"  transport  coefficients,  so  it 
seems  more  reasonable  to  use  the  simpler  Chapman-Korst  model  at  the 
present  time.  If  this  procedure  is  adopted,  then  what  remains  is  to 
develop  a  reasonable  empirical  model  for  the  eddy  coefficients.  A  major 
part  of  the  present  paper  is  devoted  to  this  task. 

Having  developed  appropriate  expressions  for  the  eddy  coefficients, 
one  can  then  solve  the  non-similar  mixing  layer  equations,  an:  thereby 
determine  how  rapidly  the  profiles  change  from  the  initial  distributions 
on  the  body  at  separation  to  the  fully-developed  profiles  far  downstream 
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from  separation.  For  the  base  flow  this  rate  of  change  or  build-up  of 

the  profiles  may  be  of  crucial  importance  if  recompression  occurs  before 

the  asymptotic  or  fully  developed  condition  is  approached.  This  is  the 

(2  3) 

situation  occurring  in  laminar  case,  *  and  could  also  occur  in  the 
turbulent  case. 

The  present  report  contains  three  separate  parts.  The  first  shows 
that  experimental  turbulent  mixing  data  can  be  correlated  by  an  eddy 
viscosity  model  based  on  the  assumption  that  pc  =  function  of  >:  only. 
Thus  the  effect  of  compressabil ity  is  taken  into  account  in  the  same 
manner  as  with  laminar  flow;  the  effect  of  turbulence  thereby  being 
placed  in  the  streamwise  variable  transformation.  The  second  part  of 
the  report  deals  with  the  solution  of  the  non-similar  mixing  layer 
problem,  assuming  that  the  previously  developed  expression  for  the 
eddy  viscosity  remains  valid  even  near  the  separation  point.  The  third 
portion  of  the  report  utilizes  the  results  of  the  non-similar  mixing 
layer  analysis  and  the  Chapman-Korst  base  flow  model  to  predict  the 
base  flow  properties  of  a  highly  cooled  slender  cone  under  re-entry 
conditions . 

As  might  be  expected,  the  results  of  this  report  depend  crucially 
on  the  expression  developed  for  the  eddy  viscosity.  Unfortunately,  the 
available  data  cover  a  very  limited  Mach  number  range  and  are  for 
essentially  adiabatic  flow,  so  that  a  considerable  extrapolation  is 
involved  in  applying  the  model  to  the  re-entry  case.  Thus  the  final 
base  flow  results  presented  must  be  considered  tentative  and  subject 
to  revision  when  new  turbulent  mixing  data  become  available. 
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2.  EDDY  VISCOSITY 

A  fundamental  assumption  of  this  report  is  that  the  compressible 

mixing  layer  flow  can  be  transformed  to  an  equivalent  incompressible 

flow.  The  transformation  is  of  the  type  used  by  Mager . ^ ^Burggraf 

(8) 

Ting  and  Libby  'and  others,  in  which  the  stream  function  and  shear 
stress  over  a  mass  element  are  assumed  invariant.  In  the  incompress¬ 
ible  plane  the  eddy  viscosity  is  assumed  to  be  given  bv  Prandtl's 
(9) 

expression,  and  is  proportional  to  the  width  oi  t,.e  mixing  zone 
and  the  velocity  difference  across  the  layer.  Ferri^^recently  has 
argued  that  when  generalized  to  compressible  f low,  Prandtl ' s  expression 
should  be  proportional  to  the  tangential  mass  flux  difference  instead 
of  the  velocity  difference.  For  mixing  layers  with  negligible  flow 
on  the  inside  both  expressions  are  identical.  A  more  serious  defect 
is  that  the  eddy  viscosity  is  assumed  to  suddenly  change  from  a  form 
suitable  to  the  body  boundary  layer  to  that  characteristic  of  the 
mixing  layer  after  the  shoulder  expansion. 


With  these  assumptions 
is  the  following: 


the  resulting  compressible  eddy  viscosity 


(1) 


In  Eq.  (1)  k  is  a  dimensionless  proportionality  factor  to  be  deter¬ 
mined  by  correlating  the  theory  with  experimental  data  and  Is  at  most 
a  function  of  the  external  flow  conditions.  The  quantity  b  is  a 
suitably  chosen  width  of  the  mixing  layer  and  is  a  reference 
density  used  to  give  e  the  dimensions  of  viscosity.  As  pointed  out 
by  Ting  and  Libby,  Eq.  (1)  implies  that  e  is  variable  across  the  mix¬ 
ing  zone  because  of  the  density  variation,  whereas  in  the  incompressible 
case  it  is  constant.  For  the  mixing  layer  the  width  is  assumed  to  be: 
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00 


u*U  .-u*) 


(2a) 


so  that  e  is  proportioral  to  the  momentum  thickness  of  the  mixing 
layer.  We  could  also  try  for  the  width 


where  (by  convention)  y.  and  y  are  respectively  the  values  of  y  where 
2  1  ^ 

u*  =  0.1  and  0.9.  It  turns  out  that  Eqs.  (2a)  and  (2b)  produce 
essentially  the  same  final  results  for  a  mixing  layer  of  negligible 
initial  thickness.  For  the  more  general  problem  however,  Eq.  (2a) 
seems  more  reasonable  than  (2b).  The  only  justification  for  these 
assumptions  for  b  at  present  is  that  they  simplify  the  analysis  some¬ 
what  by  uncoupling  the  energy  terms  from  the  momentum  equation.  As 
with  all  turbulent  theories  of  this  type,  further  justification  for 
the  selection  of  the  expression  for  e  must  come  by  comparing  the 
results  of  the  analysis  with  experimental  data. 
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3.  CONSERVATION  EQUATIONS 

The  boundary  layer  equations  involving  the  mec.n  flow  variables  are 
assumed  to  be  appropriate  for  the  description  of  the  constant  pressure 
mixing  region  (Fig.  1),  provided  the  exchange  coefficients  are  replaced 
by  their  empirical  turbulent  counterparts.  In  spite  of  the  well  recog- 
nitud  fact  that  er  rgy,  species,  and  momentum  transport  may  occur  by 
different  mechanisms  in  turbulent  flow,  we  simpl:.  Fy  matters  by  assuming 
that  Fr  =  Let=  1.  Then  only  the  eddy  visccs  .ry  is  used  explicitly. 

This  rough  model  later  may  be  refinet  the  results  prove  prot  'ing. 
The  equations  are: 


(pur  ^)  +  (pvr  ^)  =0 

o  x  o  /y 

(3a) 

puux  +  pvuy  -  (euy)y 

(3b) 

puH  +  pvH  =  (pH  ) 
x  y  y  y 

(3c) 

where  rQ  ~  rQ(x)  t^ie  mean  radius  of  the  thin  mixing  layer  in  the 
axi symmetric  case.  The  boundary  conditions  are  (Fig.  1): 


y  =  -co:  u  =  0,  H  =  Hc  (constant)  (Aa) 

y=oa:  u  =  u  ,  H  =  H  (Ab) 

e  e 

y  =  0:  v  =  0  (dividing  streamline)  (Ac) 

x  =  0,  y  >  0:  u  =  given  initial  velocity  profile  (Ad) 

H  =  H  -  (H  -  Hw)  (1-u*)  (Crocco  Integral)  (Ae) 
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Since  we  wish  to  consider  c.  es  where  Hc  t  H  ,  the  Crocco  integral 

is  not  applicable  in  the  mixing  region,  although  it  is  assumed  to  apply 

(ll) 

initially  (at  the  body  base,  x  =  0) .  Thus,  as  in  the  laminar  case, 
introduce  an  auxiliary  enthalpy  function  W  which  is  related  to  H  as 
f ol lows : 

H»H-(H-H)(1-  u*)  -  (H  -  H  )  W  (5) 

e  0  c  cw 

Then  the  W  equation  is 

put’  +  pvW  =  (rW  )  (ba) 

x  y  y  y 

W(-cs)  =  W(»)  *  0  (6b) 

W(x  =  0)  =1  -  u*  (6c) 

This  allows  the  energy  equation  to  be  solved  without  explicit  knowledge 
of  the  core  enthalpy  Hc.  With  W  available  the  enthalpy  profile  may  be 
calculated  from  Eq.  (5)  after  is  specified. 

In  order  to  transform  the  equations  to  incompressible  form,  the 
usual  transformation  of  the  y  coordinate  is  introduced.  Let 


and  define 


Then 


v  = 


pur 
e  e 


4 


u*  =  u/u 

e 

F  _  2iL* 
F  "  SY 


pe  =  k  p  r 
H  he  o 


where 


4>(x) 


o 


u*  (1-u*)  du* 
F 


(7a) 


(7b) 

(7c) 

(8a) 


(8b) 
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Introduction  of  relations  (7)  and  (8)  in  the  conservation  equations  (3) 
and  further  transformation  to  Crocco  coordinates  yields 


SF 

5x 


=  (p  u  r  Jk  <i>)  F‘ 


e  2  o 


du*2 


(9) 


and  a  similar  expression  for  W.  The  factor  in  parentheses  on  the  right 
hand  side  depends  on  x  only,  although  depends  upon  the  solution. 
Therefore  Eq,  (9)  can  be  put  in  the  form  of  the  laminar  momentum  equa¬ 
tion  if  we  define 

dS*  =  F,  (p  u  r  ^k  <I>)  dx 


w 


e  e  o 


(10a) 


F*  =  F  ''F 


w 


(10b) 


where  F^  is  the  value  of  F  on  the  body  before  separation. 
Then  the  conservation  equations  become 


dF*  2  VT* 

U'*S*  *u*2 


(Ha) 


*  3W  _  _,2  S2W 


^  =•  p* 

3S* 


(lib) 


These  equations  are  in  a  form  identical  to  the  laminar  equations  and 
also  have  the  same  boundary  conditions  as  in  the  laminar  case.  The 
differences  are  in  the  inversion  of  S*  back  to  physical  space  and  in 
the  initial  conditions.  If  the  body  boundary  layer  is  turbulent, 
initial  conditions  for  the  mixing  layer  must  take  this  into  account. 

On  the  other  hand  it  is  conceivable  that  the  body  boundary  layer  might 
be  laminar  with  transition  occurring  at  the  shoulder.  In  this  case 
the  initial  profile  wi1!  also  correspond  to  the  laminar  case. 
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For  a  ls.ainar  body  boundary  layer  the  Initial  and  boundary  condi 
tion  are; 


S*  =  0: 

F*(0,u*) 

- f VV/fY°> 

(12a) 

W(0,u*) 

=  1-u* 

(12b) 

u*  =  0: 

F*(S*,0) 

=  W(F*,0)  =  0 

(12c) 

u*  =  1 : 

F*(S*,1) 

=  W(S*,1)  =  0 

(12d) 

where  f"  is  obtained  from  the  Blasius  solution  (see  Section  7). 

b  Cl) 

These  equations  have  been  previously  solved  and  the  solution  tabulated. 

The  problems  remaining,  then,  are  the  inversion  of  the  turbulent 
flow  transformation  S*  back  to  physical  space  and  establishment  of 
initial  conditions  for  turbulent  flow  on  the  body. 
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4 .  SIMILAR  SOLUTION 


A  first  step  in  inverting  the  turbulent  solution  is  to  obtain  the 

parameter  k  from  experiments.  For  this  purpose  we  investigate  similar 

solutions  which  can  then  be  compared  with  experimental  data.  It  is 

known  from  the  laminar  case  that  at  large  values  of  S*  the  effect  of 

the  initial  structure  of  the  shear  layer  will  disappear  and  the  solution 

(  ) 

will  asymptotically  approach  the  similar  solution  of  Chapman. 

Since  the  turbulent  problem  has  been  made  mathematically  equivalent  to 
the  laminar  case  with  the  exception  of  initial  conditions,  the  asymp¬ 
totic  solutions  must  be  identical. 


In  order  to  obtain  the  Blasius  equation  instead  of  its  counterpart 
in  Crocco  coordinates  we  introduce  the  normal  distance  parameter 


n 


TF 


(2S*) 


1/2 


(13a) 


Then  the  velocity  is  assumed  to  be  a  function  of  q  only 

u*  =  f  '  Ctj  )  (13b) 

When  Eqs.  (l3a)  and  ( 13b)  are  introduced  in  the  conservation  equations 
the  familiar  results  are  obtained. 


f  +  ff  =  0 

Wi  m 


(14a) 


W  +  fW  =  0 


(14b) 


The  boundary  conditions  are  f  (-<«)  =  0,  f  (od)  =  1,  W(-oo)  =  W(x)  =  0, 
and  f(0)  =  0.  The  last  condition  follows  from  the  specification  that 
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rj  —  0  corresponds  to  the  dividing  streamline.  Clearly  the  solution  for 
W  is  W  =  0,  which  of  course  implies  the  usual  Crocco  integral  relation 
for  the  total  enthalpy  from  Eq.  (5).  The  solution  for  f  may  be  obtained 
from  tabulations  by  Chapman and  Christian  ^  \ 

The  momentum  thickness  function  for  the  similar  solution  is  given  by 

1/2 


$ 


(15a) 


w 


where 


f'(l-f')  dn(  =  ,8756 


(15b) 


Use  of  Eq,  (15)  in  Eq,  (10)  and  integration  yields; 

(2S*)^2=  k  c  F  \  p  u  r  ^dx 
w  \  e  e  o 


(16a) 


Hence 


h  = 


(16b) 


where 


cx 

\  p  u  r 


§  =  kc\  pur  Jdx  =  (2S*)^2/F 
j  e  e  o  w 

o 


(16c) 


It  can  be  seen  from  Eqs.  (16b)  and  (16c)  as  well  as  the  definition 
of  Y,  Eq.  (  7a),  that  the  similarity  parameter  q  becomes  proportional 
to  y/x  for  a  two  dimensional  incompressible  flow,  which  is  a  well 
known  behavior. 
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5.  CORRELATION  WITH  EXPERIMENT 


The  primary  source  of  correlation  is  obtained  from  experiments 
dealing  with  the  turbulent  flow  of  jets  exhausting  into  a  quiescent 
region  (Fig.  2).  The  region  most  applicable  to  the  problem  of  the  base 
flow  is  the  mixing  zone  between  the  jet  exit  and  about  five  diameters 
downstream.  This  type  of  flow  is  "inside-out"  compared  with  the  base 
flow  problem,  since  the  high  velocity  flow  is  inside  the  mixing  region 
rather  than  outside.  This  hopefully  will  not  make  any  difference,  so 
that  if  the  jet  experimental  data  can  be  correlated  with  the  theory, 
then  the  theoretical  results  can  be  applied  directly  to  the  base  flow 
problem. 

The  primary  source  of  data  for  the  eddy  viscosity  correlation  of 

(14) 

this  report  is  the  work  of  Maydew  and  Reed,  who  measured  velocity 

profiles  in  the  near  mixing  region  of  turbulent  jets  exhausting  into 

the  atmosphere.  The  nozzle  exit  diameter  of  these  experiments  was  3", 

and  profile  data  were  obtained  for  exit  Mach  numbers  .7,  .85,  .95,  1.49, 

and  1.96  at  five  axial  stations  in  the  range  1,5"  <  x^ll.5".  Leipmann 

(15) 

and  Laufer  also  give  data  for  Mfc.05,  while  some  earlier  data  of 
possible  applicability  to  the  present  problem  are  referred  to  by  Maydew 
and  Reed, 

One  result  which  can  be  compared  with  the  theory  is  the  velocity 

on  the  dividing  streamline.  According  to  the  solution  of  Eq,  (14«)  this 

value  should  be  u*  =  .587  if  the  flow  is  sufficiently  downstream  of 

separation  for  similarity  to  hold.  In  order  to  check  this  result,  the 

(14) 

jet  mixing  profiles  of  Maydw  and  Reed  were  integrated  outward  from 
the  axis  of  the  jet  until  the  mass  in  the  profile  was  equal  to  the  mass 
flowing  throe  Ne  jet.  In  these  calculations  the  Crocco  integral 
relation  was  used  to  calculate  the  density  variation,  and  the  entrance 
mass  was  calculated  from  the  measured  stagnation  temperature  and 
pressure,  assuming  an  isentropic  one-dimensional  expansion.  Further 
details  regarding  this  mass  balance  are  given  in  the  Appendix. 
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Fig.  3  shows  the  calculated  results  of  the  mass  balance  for  the  five 
exit  Mach  numbers  of  the  experiments  of  Maydew  and  Reed.  One  sees  that 
the  calculated  experimental  dividing  streamline  velocity  u*  is  about 
0,6  or  perhaps  a  little  higher.  It  is  substantially  independent  of 
streamwise  distance  for  x/d  ^  1,  indicatiig  that  the  mixing  region  has 
probably  reached  similarity.  It  should  be  noted,  however,  that  the 
accuracy  of  the  integration  procedure  is  very  poor  for  stations  close 
to  the  jet  exit,  since  most  of  the  mass  flow  inside  the  dividing  stream¬ 
line  is  in  the  potential  core  and  not  in  the  mixing  profile  (see  Fig.  2). 
Under  these  conditions  a  small  error  in  the  mass  flow  calculation  will 
produce  a  large  error  in  the  value  of  u*  on  the  dividing  streamline. 
Further,  the  use  of  the  Crocco  integral  for  the  density  variation  may 
be  questioned.  We  therefore  conclude  from  Fig.  3  that  the  dividing 
streamline  velocity  data  agree  with  the  present  theory  within  the 
possible  error  of  the  measurements  and  integration  procedure,  but  these 
data  are  probably  too  crude  to  constitute  a  critical  test. 

The  second  and  more  important  portion  of  the  data  correlation  concerns 
the  fitting  of  the  theoretical  velocity  profile  to  the  data.  This  has 
traditionally  been  done  by  finding  the  best  numerical  value  of  the 
"spreading  parameter"  o  such  that  the  velocity  profile  is  represented  by 

U*  =  g(c  h  (17) 

where  g  is  a  profile  function  specified  beforehand.  For  example,  the 
error  function  profile  is  frequently  used,  while  Maydew  and  Reed  find 
that  the  data  correlate  well  with  the  results  of  Crane. The  numeri¬ 
cal  value  of  o  for  "best  fit"  of  course  depends  on  the  choice  of  profile. 

A  certain  arbitrariness  also  exists  in  the  selection  of  the  profile 
position  from  which  y  is  measured,  and  this  traditionally  is  taken  as 
the  location  where  u*  =  . 5  for  the  data. 
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We  begin  by  noting  that  the  present  solution  for  compressible  flow 
does  not  allow  the  physical  velocity  profile  to  be  expressed  in  the  form 
of  Eq.  (17).  This  of  course  is  due  to  the  manner  in  which  the  effect 
of  compressibility  is  included  in  the  transformation.  To  derive  the 
form  of  the  velocity  profile  in  the  physical  coordinates,  the  Howarth 
integral  in  the  expression  for  "H  must  be  inverted.  This  can  be  done 
by  assuming  perfect  gas  behavior  (which  is  reasonable  for  the  Maydew- 
Reed  experiments)  and  using  the  Crocco  integral  relation.  Thus  for  a 
perfect  gas; 


h_ 

h 


K 


e  L_  e 


(18a) 


H  .  9 

-—  =  1  +  ^  M 
h  2  e 

c 


(18b) 


Then  if  r  is  assumed  to  be  constant,  one  finds  that 
o 


<  rt 


1  - 


H 


2 

f  dn 


(19) 


Since  u*  =  f^O")),  it  is  clear  that  an  expression  for  u*  of 
Eq.  (17)  is  not  possible. 


the  form  of 


To  compare  with  the  velocity  profile  data  of  Maydew  and  Reed, 

Eq.  (19)  is  multiplied  by  a,  where  a  is  now  considered  to  be  just  a 
given  scale  factor  for  each  set  of  profile  data.  Also  In  the  above, 
for  convenience  y  is  measured  from  the  position  for  u*  =  .587,  rather 
than  u*  =  .5,  as  done  by  Maydew  and  Reed. 

It  is  assumed  that  H^/H^  =  1,  since  this  was  assumed  in  the  Maydew- 

Reed  data  reduction,  even  though  .93  <  H  /H  <  1  for  their  experiments. 

c  e 
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Also  x  =  x  +  Ax,  where  x  Is  the  experimental  distance  from  the  jet 
e  e 

exit  and  Ax  is  an  additional  incremental  d’st.  nee  to  the  "virtual  origin" 

(14) 

of  the  turbulent  mixing  layer  as  tabulated  by  Maydew  and  Reed, 

Fig.  4  shows  the  correlation  of  the  Maydew-Reed  data  with  the  Chapman 
profile  for  the  various  experimental  Mach  numbers.  The  constant  k  was 
first  obtained  by  finding  the  "best  fit"  of  Eq,  (19)  with  each  set  of 
profile  data,  using  a  "le>  t  squares"  analysis,  (Note  that  the  momentum 
thickness  factor  c  =  ,8756  from  the  Chapman  solution.)  The  theoretical 
curves  shown,  however,  are  based  on  the  final  equation  for  the  correla¬ 
tion  constant,  as  described  below.  It  is  seen  that  the  shapes  of  the 
profiles  seem  to  correlate  quite  well.  This  is  perhaps  not  too  surpris¬ 
ing,  however,  since  many  curves  of  this  general  shape  seem  to  correlate 
in  a  reasonable  manner  if  the  adjustable  constant  is  properly  chosen. 

It  was  hoped  initially  that  the  quantity  k  would  not  be  a  function 
of  Mach  number,  but  would  be  a  universal  constant.  Instead  it  is  found 
to  vary  quite  strongly  with  Mach  number  for  .05  5  M  s  3,  It  is  reason¬ 
able  to  assume  that  this  variation  is  related  to  the  density  ratio 
across  the  mixing  region  p^/p^.  For  a  perfect  gas  this  is  equal  to 

H  /h  ,  and  if  the  flow  is  completely  adiabatic,  then  p  /p  =  H  /h  = 
c  6  .  a  e  c  e  e 

y  —  1  2. 

(1  +  Mg  ).  Figure  5  shows  a  log-log  plot  of  k  as  a  function  of 

Pe/Pc,  assuming  adiabatic  flow.  Based  only  on  the  Maydew-Reed  data, 
it  is  found  that  very  nearly 


If  the  data  of  Leipmann  and  I.aufer  and  tentative  data  of  Zumwalt 
(see  Appendix)  had  been  included  in  the  correlation,  the  constants  in 
Eq,  (20)  would  probably  be  changed  slightly.  Thus  for  the  similar 
solution:  /  v  b. 


C  C  =  b 


/  \ 

2  (°c\  2 

lpe  \r)  u< 
\  e/ 


(21) 


Comparing  this  with  Eq,  (1),  it  is  seen  that  a  more  reasonable  reference 

2  2 

density  probably  would  have  been  rather  than  . 
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6.  NGN-SIMILAR  SOLUTION 

Equation  (21)  is  a  plausible  form  for  the  eddy  viscosity,  which 
should  be  applicable  at  least  to  fully-developed,  .‘.early  adiabatic 
turbulent  mixing  for  Mach  numbers  up  to  about  3.  The  justification  for 
its  use  lies  in  the  fact  that  experimental  data  correlate  well  with  theory. 
Application  of  this  equation  at  higher  Mach  numbers,  for  highly  cooled 
mixing  layers,  or  in  the  non-similar  region  close  to  the  separation  point 
(Fig.  1)  obviously  may  be  impropei .  Because  of  the  urgent  need  for  some 
sort  of  theoretical  description  of  turbulent  mixing  under  these  condi¬ 
tions,  however,  the  extended  use  of  Eq.  (21)  is  proposed.  Any  conclusions 
based  on  this  model  must  of  course  remain  tentative  for  conditions  out¬ 
side  of  the  established  correlation  range. 


We  now  propose  to  solve  non-similar  Eqs.  (  11)  using  Eq .  (2l)  for  the 
eddy  viscosity.  We  introduce  the  correlation  in  the  x  coordinate  trans¬ 
formation  for  the  non-similar  case. 


S* 


w 


peWo' 


dx 


where 


bo 


PT  = 


— 

fo  \  1  1 

1 

—  i  $ 

1 

=  k  $> 


(22a) 


(22b) 


(  3^ 

Comparing  Eq.  (22^  with  the  well-known  laminar  expression  for  S*,  ' 

we  see  that  the  term  replaces  the  laminar  viscosity  ti  times  r 

i  e  o 

Thus  the  relative  growth  rates  dS-'/dx  of  turbulent  and  laminar  mixing 
layers  (in  the  two-dimensional  case;  are  directly  related  by  the  ratio 
.  One  sees  that  if  the  eddy  viscosity  is  in  fact  representable 
by  Eq.  (21)  over  a  wide  ran^e  of  conditions,  then  the  density  ratio 
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across  the  mixing  layer  (p  / p  )  is  extremely  important  in  determining  the 
growth  of  the  turbulent  mixing  ^ayer.  It  has  already  been  shown  that 
the  conservation  equations  have  been  converted  to  laminar  form  in  Eqs.  (ll). 
Eqs.  (22)  provide  the  transformation  of  S,c  back  to  physical  space.  In 
essence  what  has  been  accomplished,  therefore,  is  the  placement  of  the 
turbulent  effects  into  the  coordinate  transformations.  One  should  thus 
be  able  to  obtain  universal  solutions  to  Eq.  (11)  which  depend  only  on 
the  shape  of  the  initial  profile  F(0,u*),  just  as  in  the  laminar  case. 

The  task  that  remains  is  to  obtain  the  initial  mixing  layer  profiles 

by  solving  the  boundary  layer  on  the  body  surface  upstream  of  separa- 

(2  3) 

tion.  Eqs.  (11)  will  then  be  solved  by  a  finite  difference  method.  ’  ’ 

The  simplest  case  to  consider  is  that  for  laminar  flow  on  the  body  sur¬ 
face,  which  is  a  limiting  situation  which  might  correspond  to  the  occur¬ 
rence  of  transition  at  the  separation  point.  A  slightly  more  difficult 
but  still  tractable  problem  occurs  for  fully  developed  constant  pressure 
turbulent  flow  (i.e.  cone  or  wedge).  The  calculations  of  this  report 
are  restricted  to  these  cases. 
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7.  LAMINAR  INITIAL  CONDITIONS 

For  laminar  flow  on  the  body  it  is  found  that  Elasius  solution  is 
applicable.  Thus  ac  the  separation  point: 


F(0,u*)  = 


(23a) 


I* 


-  f  VV 


u  r  1  C 
e  o  \ 

\[2S  ° 

v  w 


pdy 


(23b) 


(23c) 


21 


S  =  \  C,  p  u  u  r  Jdx 

w  \  b  e  e  e  o 


body 


(23d) 


C,  =  pp/p  p  (Chapman-Rubesin  constant 
6  G  for  the  body) 


(23e) 


In  (23a)  and  0.3b)  prime  denotes  differentiation  with  respect  to 

(2) 

the  Elasius  variable  n  .  As  in  the  laminai  mixing  case  we  let: 

D 


F  =  f"  (0)/\l  2S 
w  B  '  v  w 


(24a) 


ff,B(0)  =  .4696 


(24b) 


The  initial  condition  for  the  energy  equation  is  W  =  1-u*,  since 
the  Crocco  integral  is  assumed  to  apply  to  the  flow  on  the  body  surface. 
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8.  TRANSFORMATION  INVERSION  -  LAMINAR  BODY 

The  calculation  of  S*  is  now  considered.  From  Eq,  (22),  we  are  led 
to  the  following  differential  equation: 


o 


Define  a  new  variable  x*  by  the  following  relation 


x* 


r 

\  p  u  F  k  r  "*dx 
J  e  e  w  o 


(25a) 


(25b) 


(26) 


Then  Eq.  ( 25a)  may  be  written 


dS* 

dx* 


4*v(s*) 


S*(0)  =  0 


(27) 


This  equation  was  integrated  using  the  function  data  for  from  the 

shear  layer  solution,  and  the  resulting  function  x*(S*)  is  given  in 
Fig.  14.  The  relation  between  S*  and  physical  x  is  therefore  obtained 
from  Eq.  (26)  and  the  function  xy'(S*)  . 

For  comparison  of  turbulent  mixing  with  the  laminar  mixing  case, 
it  may  be  desirable  to  express  the  turbulent  solution  in  the  laminar 
variables.  For  this  we  need  to  calculate  S**  for  the  turbulent  case, 
where 

S**  =  S/S 

w 


S 


f 

Jo 


2 4 

C  p  u  u  r  Jdx 
m  e  e  ®  o 


(28) 
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and  C  is  the  laminar  Chapman-Rubes  in  constant,  for  the  mixing  laver. 
m 

Then  the  rate  of  build-up  of  transformed  length  scale  ratio  S**  with 
respect  to  x*  is  given  by: 


dS** 

dx* 


k  F  S 


w  w 


(29) 


The  bracket  term  in  Eq.  (29)  is  a  constant  for  a  given  body  at  given 
flow  conditions.  The  value  of  this  constant  in  a  base  flow  depends  on 
the  details  of  the  "matching,"  which  is  carried  out  using  the  Chapman - 
Korst  recompression  condition.  For  two  dimensional  flow  the  laminar 
length  variable  S**  is  proportional  to  x*(S*),  but  for  axisymmetric 
flow  the  radius  factor  enters  explicitely.  Further  discussion  of  the 
laminar -body ,  turbulent-mixing-layer  problem  is  deferred  until  later. 
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9.  TUPBULENT  INITIAL  CONDITIONS 


For  a  turbulent  body  boundary  layer  we  need  to  obtain  new  initial 
conditions  which  replace  Eq.  (23)  in  the  analysis  of  the  preceding 
sections.  In  spite  of  the  great  effort  which  has  been  devoted  to 
research  on  turbulent  boundary  layers,  there  still  does  not  exist  a 
well  established  theoretical  method  to  obtain  compressible  turbulent 
boundary  layer  profiles.  The  most  recent  and  promising  method  appears 
to  be  the  transformation  method  of  Coles  ,  ^  1  as  further  explained  by 

(  I  QJ 

Crocco.  '  Thus  a  transformation  will  be  found  which  (hopefully) 
establishes  a  correspondence  between  a  known  incompressible  flow  and 
the  desired  compressible  flow.  The  incompressible  profiles  (which  are 
established  by  a  semi -empirical  method)  are  thereby  transformed  to  the 
compressible  flow. 

Let  "barred"  symbols  refer  to  the  transformed  incompressible  flow 
and  unbarred  symbols  represent  the  compressible  flow.  The  aim  of  the 
Coles  transformation  is  to  find  the  quantities  c(x)  ,  r|  (x) ,  and  §(x) 
defined  as  follows : 


^00 


(30a) 


_i_  = 

p  dy 


(30b) 


(30c) 
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vheie  the  stream  functions  are: 


<30d) 


j  fy 

Jr  =  r  .  pu  dy 

°  J 


(30e) 


(The  quantity  of  cr  Eq.  (30a)  is  not  related  to  the  jet-mixing  c  of 
Fig.  4.)  Restricting  attention  to  the  case  of  zero  pressure  gradient. 
Coles  and  Crocco  find  that 


u 

u 


constant 


The  relation  between  the  wall  shear  stresses  is 


on 


T  p  U  C-P  U 

w  c  v  w  _  fee 

x~  2pu_„  2 

W  T]  W  W  urp  U 

fee 


(32) 


? 

where  Cc=  x  / (p  u  ~/2) 
f  wee 

introduces  the  idea  of 
the  result 


is  the  usual  skin  friction  coefficient.  Coles 
a  turbulent  substructure,  which  yields  for  c 


(33) 


where  ug  is  a  mean  substructure  viscosity,  obtained  from  the  hypothesis 
of  a  constant  substructure  Reynolds  number .  ’  ^Vrom  Eq,  (3l)  to  (33) 
we  obtain 


2 


(34) 
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For  the  case  of  Pr  =  1  the  Crocco  integral  relation  holds  on  tin  body. 

Then  if  laminar  viscosity  is  assumed  proportional  to  static  enthalpy, 

the  following  generalization  of  Coles'  results  (Coles'  Fq.  4.17  for 

u  /u  is  obtained: 
s  w 


a  «  17.5 


a2=  305 


(35b) 


The  simplest  way  to  treat  the  compressible  turbulent  boundary  layer 
using  this  formulation  is  to  specify  C^  (i.e.  work  the  problem  backwards 
by  specifying  the  equivalent  incompressible  skin  friction  coefficient). 
Then  the  compressible  skin  friction  is  given  by  Eq.  (37),  and  the  com¬ 
pressible  heat  transfer  is  obtained  from  the  enthalpy  gradient  using 
the  Crocco  integral.  What  remains  to  be  done,  then,  is  to  find  the 
relation  between  the  compressible  and  incompressible  length  scales 
(i.e.  Reynolds  numbers).  Then  the  compressible  velocity  profile  can 
be  expressed  in  terms  of  a  given  incompressible  profile  corresponding 
to  the  specified  C^. 

Rather  than  finding  the  ’•elation  for  §(x),  we  instead  seek  a  direct 
relation  between  C^  and  the  compressible  length  x.  If  the  Coles-Crocco 
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transformation  relations  are  substituted  in  the  compressible  momentum 
integral  relation,  it  is  found  that 


-’.x 


pur  JJx 
w  e  o 


R 

e 


e 


u*)  dy 


(38a) 


(38b) 


The  incompressible  momentum  thickness  Reynolds  number  R  is 

e9 

related  to  by  a  relation  based  on  the  law  of  the  wall.  Let  the  in¬ 
compressible  velocity  profile  in  this  region  be  given  by  (Ref.  19, 
page  140) : 


u_ 

u 


(39a) 


u 


(39b) 


where  according  to  Coles, 2<^K  =  .4  and  C  =  5.1.  Then  from  Eq.  (38b) 
it  is  found  that 


-KC 

e _ 

K 


(1  -  T>  +  (1 


(40a) 


where 


(40b) 


By  using  Eq.  (40)  and  (35),  Eq.  (38)  can  be  integrated  along  the  body 
from  the  point  of  transition  to  the  uase.  If  £  is  assumed  to  be  very 
large,  then  the  integration  is  easy  to  carry  out.  Assuming  that  transi¬ 
tion  occurs  at  the  nose,  the  final  result  may  be  put  into  the  standard  form 
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=  ri  +  C2  lGglO 


Vw 


C,  R 


r  L\ 

u  o  f  e 
e  s  T 


where  cp  is  defined  in  Eq.  (34)  and 
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(Ala) 


(alb) 


(41c) 


(Aid) 


Equations  (35),  (37),  and  (41)  are  three  equations  for  the  three 

—  dy 

unknowns  C,f  C,  ,  — .  The  parameters  are  H  /H  ,  H  =h  ,  and  R  ,  which 

f  f  4  s  we’ee  eT 

will  be  known  for  a  given  body  at  given  flight  conditions. 


For  wedges  and  cones  Eq,  (4lc0  may  be  integrated  explicitely. 
Figs.  6,  7  and  8  show  values  of  the  equivalent  for  these  bodies  as 
a  function  of  wall  enthalpy  ratio,  boundary  layer  edge  Mach  number 
or  enthalpy  ratio  H^/h^,  and  reduced  Reynolds  number  R: 


R 


C  Re. 
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From  the  result  for  C  ,  the  initial  compressible  F(u,v)  profile  at  the 
body  base  can  be  found  from  the  laminar  sublayer  and  law-of -the -wall 
prof i les ,  Thus  : 


Sublayer :  F 


i  r  J 

w  o 


(43a) 


in  the  region 


0 


i* 


+ 


Nf 


Cf/2 


(43b) 


Rubesin^^f  inds  that  u+=  13.1,  In  the  law-of -the-wall  region  the 
velocity  gradient  function  is  found  from  Eq.  (39)  and  the  proper*  -‘S 
of  the  Coles-Crocco  transformation: 


Law-of -the-wall : 


(43c) 


in  the  region 

u+  \|"£f/2  <  u*  <  1 


(43d) 


In  this  form  the  initial  velocity  gradient  profile  l/F^  is  a  func¬ 
tion  only  of  the  equivalent  incompressible  skin  friction  coefficient  C^, 
Therefore  this  is  the  natural  parameter  defining  a  family  of  solutions 
to  the  non-similar  turbulent  mixing  layer  problem.  Fig.  9  shows  the 
shape  of  these  profiles  for  various  values  of  in  the  range  of  interest. 
One  sees  that  there  are  square  corners  at  locations  corresponding  to  the 
sublayer  limit,  and  the  outer  edge.  This  behavior  is,  of  course,  not 
physically  reasonable,  so  these  points  should  be  arbitrarily  rounded 
slightly  to  obtain  a  reasonably  smooth  starting  profile  for  the  non¬ 
similar  mixing  layer  calculation.  This  hopefully  will  not  affect  the 
Dividing  streamline  properties  much  because  the  dividing  streamline  is 
initially  at  u,v  =  0,  One  could  probably  improve  the  profile  by  includ¬ 
ing  the  law-of -the-wake  and  buffer  layers,  but  this  does  not  seem  worth 
the  effort  at  presem. .  The  discontinuity  at  S*  =  u*  =  0  (i.e,  F(0)  =  0 

immediately  after  separation)  is  of  course  treated  in  the  same  manner 

(22) 

as  done  in  the  laminar  case.  s 
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Note  that  for  a  base  flow  p-oblem  the  sudden  turn  at  the  separation 
point  causes  a  distortion  of  the  profile.  If  one  assumes  that  this  occurs 
according  to  an  isentropic  expansion  along  streamlines  (as  in  Ref.  23), 
then  one  could  use  the  distorted  initial  profiles  as  the  mixing  layer 
initial  conditions.  This  would  involve  finding  the  streamlines  using 
the  Coles-Crocco  relation  for  the  stream  function  from  Eq.  (30)  and 
presents  no  essential  difficulty.  The  presumed  increase  in  accuracy 
from  this  refinement  does  not  appear  to  be  worth  the  effort  at  this 
stage,  but  can  be  included  in  a  later  analysis  if  the  present  results 
appear  promising. 

For  the  turbulent  body  and  tu.’ou’ent-mixing-layer  problem,  then,  the 
initial  conditions  are  given  by: 


r 


F*  = 


0  <  u*  <  13,1  \  Cf/2 


(44a) 


exp 


r—  — 

-.4 

u* 

-  5.1 

1 

04 

(4-4 

tu 

_l 

13.1^  Cf/2  <  u*  <  1 


(44b) 


The  integration  procedure  is  again  identical  to  the  purely  laminar  case. 
The  non-similar  calculation  of  the  shear  layer  must  of  course  be 
repeated  because  the  turbulent  initial  profile  shapes  are  different. 
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10.  TRANSFORMATION  INVERSION  -  TURBULENT  BODY 

As  before,  one  needs  to  know  the  relation  between  S*  and  the  physical 
length  x.  Also  for  comparison  with  the  purely  laminar  results  it  may  be 
desirable  to  know  the  solution  in  terms  of  the  laminar  streamwise  scale 
parameter  S**  =  S/S^.  Examination  of  the  transformation  equations  reveals 
that  the  rate  of  build-up  of  S**  along  the  shear  layer  is  given  by  Eq.  (29), 
with  x*(S*)  defined  by  Eq.  (27).  Thus  the  difference  between  turbulent 
mixing  cases  having  laminar  or  turbulent  initial  conditions  resides  only 
in  the  shape  of  the  initial  profile  and  in  the  value  of  wall  velocity 
gradient  at  separation  F  .  For  the  laminar  body  this  quantity  is  given 
by  Eq.  (24),  while  for  the  turbulent  body  it  is  obtaine  '  from  Eq.  (4ia). 


EOS  RN-24 ,  5-65 


27 


li.  NUMERICAL  SHEAR  LAYER  CALCULATIONS 


Non-similar  shear  layer  calculations  were  carried  out  corresponding 

to  the  three  turbulent  initial  profiles  and  the  laminar  profile  shown  ir 

Fig.  Th.  calculations  were  started  at  Sv"  *  10  ^  and  carried  to 

3 

S*  =  10  ,  with  output  at  intervals  of  0.2  in  log^CS*).  Approximately 
ten  minutes  of  IBM  7094  time  was  required  for  each  initial  profile. 

The  u*  mesh  contained  80  intervals  in  the  range  0  s  u*  s  .05  and  152 
intervals  in  the  remaining  rang'  .05s  u*  s  1.0,  Overall  integral 
balances  (from  momentum  and  energy)  agreed  within  0. 5°/o  for  all  condi¬ 
tions,  which  is  probably  indicative  of  the  accuracy  of  the  numerical 
calculations. 

Fig,  10  shows  the  results  of  the  non-similar  shear  layer  calcula¬ 
tions  for  the  velocity  gradient  function  F*  for  the  four  different 
initial  profiles  shown  in  Fig.  9.  This  is  the  solution  of  Eq.  (11a). 

At  small  values  of  the  strearawise  variable  S*  the  profiles  resemble  the 
initial  conditions,  but  as  S*  increases  the  profiles  become  more  rounded 
and  decay  in  amplitude.  This  is  of  course  to  be  expected  since  the 
differential  equation  is  parabolic.  One  would  expect  that  as  S*  "  08 
the  shape  of  th  F*  curves  would  approach  the  asy jptotic  shape  given 

by  the  Chapman  profile.  Because  of  computer  cost,  however,  it  was 

3 

necessary  to  stop  the  calculation  at  S*  =  10  .  Only  the  Blasius  (and 
perhaos  the  C^  =  ,006)  profile  were  near  the  asymptotic  shape  at  thi3 
value  of  S*. 

The  solution  for  the  enthalpy  function  W  from  Eq.  (lib)  is  shown 
in  Fig.  11.  Bec.uuse  tne  Crocco  integral  for  total  enthalpy  is  assumed 
to  be  valid  initially,  W  =  1  -  u*  for  all  profiles  at  the  initial  sta¬ 
tion.  Since  the  W  equation  is  also  parabolic,  the  decay  of  this  function 
is  qualitatively  similar  to  the  F*  curves.  Differences  in  W  results  for 
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the  various  initial  F*  profiles  do  not  become  apparent,  however,  until 

-2 

some  distance  downstream  (S*  ^  10  )  because  all  W  profiles  begin  with 

the  same  initial  condition.  At  large  J  one  would  again  expect  the  shape 

of  the  W  function  curves  to  be  independent  of  the  initial  Fr  profile 

3 

shape,  but  this  occurs  at  S*  >  10  . 

Fig,  12  shows  the  results  for  the  velocity  u*  on  the  dividing  stream- 

°  (2) 

line,  as  obtained  by  integrating  the  momentum  equation  fot  v  =  0. 


u*  IS  . 

D  dS*  D\^u* j 


D 


(45) 


For  the  Blasius  and  C  =  .006  initial  profiles  u*  has  effectively  reached 

f  3  D 

the  Chapman  limiting  value  of  .587  at  S*  =  10  ,  but  the  other  profiles 

apparently  require  several  more  decades  in  S*  to  reach  the  limit.  All 

the  initial  profiles  used  give  the  same  value  of  u*  for  small  S*,  since 

(22) 

the  starting  profiles  of  Baum  (snail  S*  and  u*)  are  identical,  as 

may  be  seen  irom  Fig.  9. 

The  development  with  distance  of  the  enthalpy  function  on  the 

dividing  streamline  is  illustrated  in  Fig.  13.  The  limiting  value  W  =.611 

(22)  ^ 

as  S*  -  0  was  obtained  from  the  starting  solution  of  Baum.  ^  At  large 

3 

S*  the  W  function  decays  uniformly  to  zero  (Fig.  11),  but  at  S*  =  10  , 
still  has  an  appreciable  magnitude  for  the  =  ,004  and  .006  initial 
profiles. 

It  should  be  noted  that  a  direct  comparison  between  the  laminar  and 
turbulent  cases  cannot  be  made  on  the  brsis  of  Fig.  10-13,  sine,  the 
streamwise  variable  S*  is  not  directly  related  to  the  streamvise  distance. 
For  this  or.e  must  use  the  variables  x*  or  S**  for  the  turbulent  shear 
layer,  or  better  yet  the  actual  distance  x.  The  relation  between  x*  and 
3*  from  Eq.  (27)  is  shown  in  Fig.  14.  Fig.  15  gives  values  of  the  momentum 
thickness  integral  $*(S*)  which  appears  in  Eq.  (27), 
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12.  BASE  FLOW  OF  CONE  OR  WEDGE 


Fig.  lb  illustrates  the  application  of  the  mixing  layer  analysis  to 
the  base  flow  region.  The  "core"  or  recirculating  region  is  assumed  to 
have  negligible  velocity  and  constant  (but  initially  unknown)  enthalpy 
.  The  recompression  region  is  assumed  to  be  small,  and  recompression 
is  assumed  to  be  isentropic  along  streamlines.  The  distortion  of  the 
initial  profiles  at  separation  is  neglected,  although  it  could  be  included 
later  using  the  method  of  Ref.  23. 

(2  3) 

Following  the  previous  laminar  analysis,  ’  the  configuration  of  the 
base  flow  is  determined  using  the  empirical  Chapman-Korst  recompression 
condition.  This  states  that  the  total  pressure  on  the  stagnating  stream¬ 
line  just  before  reccmpression  must  equal  the  static  pressure  after  re¬ 
compression  (determined  from  the  invisci  1  flow  calculation).  Assuming 
that  the  values  of  u*  and  W  on  the  dividing  streamline  are  available 
(from  Figs.  12  and  13),  matching  involves  the  simultaneous  calculation  of 
the  inviscid  flow  (as  defined  by  the  initial  wake  -ogle),  the  core  enthalpy 
Hc ,  the  value  of  S*  ci  esponding  to  the  position  of  recompression,  and 
possibly  the  total  base  heat  transfer  rate  .  As  shown  in  Refs.  3  and  11 
these  latter  quantities  are  related  by  an  overall  energy  balance  condi- 
t  ion . 

By  equating  the  energy  entering  the  base  region  through  the  body 
boundary  layer  to  that  leaving  through  the  neck  and  by  base  heat  trans¬ 
fer,  the  following  equation  is  obtained  :  ^ 


H  =  H 
c  e 


H  Q*  +  (H  -  H  )  (K*  -  J*) 
e  b _ e  w  _ _ 

L* 


(46a) 


where 


0  F 

w 

H  (2n)J 

e 


(46b) 
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=  ♦*(()) 


(46c) 


K*  - 


-n 

\ 

J 

■  o 


du» 


body 


r-n  1 


J*  = 


u*W  .  . 
-  au“ 


recompression 


L*  = 


r-rl 


U” 

s 


u*(l-a*-W) 


F* 


recompression 


(46d) 


(46e) 


and  u*  is  the  velocity  on  the  stagnating  streamline  just  before  recom¬ 
pression.  It  may  be  seen  that  Eq„  (46a)  is  a  relation  between  Hc ,  > 
and  S*,  the  value  ot  S*  defining  the  position  of  recompression. 

Because  the  base  heat  transfer  enters  into  the  matching  analysis, 

one  needs  an  additional  relation  which  specifies  this  quantity.  Probably 

is  proportional  to  the  enthalpy  difference  (H^-  H^)  ,  where  is  the 

enthalpy  corresponding  to  the  base  wall  temperature.  An  analysis  of 
(24) 

the  laminar  case  showed  that  conditions  were  such  that  could 
safely  be  assumed  to  be  zero  without  affecting  the  base  flow  solution 
much.  For  the  time  being  it  is  assumed  that  this  conclusion  is  valid 
for  turbulent  flow  as  well,  so  that  in  what  follows  the  assumption 
=  0  is  made. 

In  order  to  compare  the  turbulent  base  flow  results  with  laminar 
results,  matching  calculations  were  carried  out  for  a  10°  cone  for  per¬ 
fect  gas  conditions  with  y  =  1.4  and  viscosity  proportion  to  tempera¬ 
ture  to  the  .76  power.  For  simplicity  the  outer  edge  conditions  in  the 
base  region  were  obtained  from  a  Frand t 1 -Meyer  expansion  at  the  corner, 
and  the  shear  layer  was  assumed  to  be  straight  between  separation  and 
recompression.  The  recompression  was  assumed  to  be  isentropic  and  all 
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Chapman-Rubes  in  constants  were  assumed  to  be  equal 
corresponds  to  calculations  previously  carried  out 
Details  of  the  numerical  procedure  for  carrying  out 
were  given  in  Ref.  11. 


to  unity.  This  model 


for  the  laminar  case, 
matching  calculations 


(25) 
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13.  RESULTS  FOR  10°  CONE 

Even  with  the  approximations  listed  in  the  previous  section,  the 
base  flow  results  are  a  function  of  Mach  number,  Reynolds  number,  body 
shape,  wall  enthalpv  ratio,  and  condition  of  the  boundary  layer  at  the 
separation  point.  Because  of  the  preliminary  nature  of  the  present 
theory  for  turbulent  flow,  it  does  not  seem  appropriate  to  make  an 
exhaustive  study  of  the  effects  of  each  variable  at  this  time.  Instead, 
attention  is  restricted  mainly  to  a  10°  half-angle  rone  with  a  "highly- 
cooled"  surface.  The  results  therefore  will  indicate  how  the  theory, 
developed  from  a  correlation  of  experimental  data  at  low  supersonic 
speeds,  is  extrapolated  to  conditions  tvpical  of  re-entry  conditions 
on  a  slender  body. 

Figure  17  shows  the  effect  of  Mach  num'  er  on  the  dividing  stream¬ 
line  velocity  just  before  recompression  for  a  cold-wall  10°  cone  at 
Re^=  10^,  The  curve  labeled  "turbulent"  corresponds  to  a  fully  developed 
turbulent  flow  on  the  body  surface  and  a  turbulent  mixing  layer.  The 

laminar  curve  gives  the  results  previously  reported  for  completely 
(25) 

laminar  flow,  while  the  "laminar-turbulent"  curve  is  based  on  the 
assumption  of  a  laminar  body  and  turbulent  mixing  layer  (i.e.  transi¬ 
tion  at  the  separation  point).  One  sees  that  at  high  Mach  numbers  all 
curves  give  a  dividing  streamline  velocity  in  the  range  0.2  to  0.3, 
which  is  very  much  lower  than  the  Chapman  value  of  .587.  This  indicates 
that  the  mixing  layer  in  the  base  flow  undergoes  recompression  long 
before  tne  fully  developed  or  asymptotic  condition  is  approached.  At 
low  Mach  numbers,  however,  the  turbulent  u*  curves  approach  u*  =  .587, 
indicating  that  at  lower  Mach  numbers  the  turbulent  mixing  effect  (as 
measured  by  the  eddy  viscosity  <e)  is  very  much  stronger.  The  strong 
Mach  number  effect  on  e  is  directly  related  to  the  effect  of  Mach  number 
on  the  eddy  viscosity  factor  k  of  Fig.  5.  In  fact,  the  effect  of  Mach 
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number  01  k  completely  overshadows  the  effect  of  Mach  on  the  initial 
profile  shape.  Thus  from  Fig.  6  it  is  seen  that  reducing  the  Mach 
number  reduces  while  Fig.  12  shows  that  this  reduces  the  rate  of 
build-up  of  dividing  streamline  velocity  u*.  Because  the  expression 
for  k  completely  dominates  the  analysis,  experiments  are  needed  to  see 
if  the  trend  of  Fig.  5  persists  to  high  Mach  number  conditions. 


Figure  18  shows  the  results  of  the  energy  balance  to  determine 
Lhe  enthalpy  in  the  recirculating  core  region.  Again  the  effect  of 
Mach  number  is  evident,  the  results  indicating  that  at  low  Mach  numbers 
the  turbulent  mixing  is  so  rapid  that  H^-*  ,  i.e.  adiabatic  conditions. 

This  occurs  in  spite  of  the  fact  that  the  upstream  body  surface  is 
highly  cooled.  Any  base  heat  transfer  would  of  course  tend  to  lower 
Hc ,  and  this  may  be  an  important  effect  in  turbulent  flow.  If  the  pro¬ 
posed  correlation  for  k  of  Fig.  5  is  correct,  then  the  core  enthalpy 
may  be  extremely  important,  since  for  a  perfect  gas; 


k 


,  06  Ou 


(47) 


The  core  enthalpy  therefore  enters  into  the  determination  of  the  eddy 
viscosity. 

The  calculated  effects  of  Mach  and  Reynolds  numbers  on  the  center- 
line  stagnation  enthalpy  at  recompression  are  shown  in  Figs.  19  and  20. 

The  results  for  h  /H  as  a  function  of  Mach  number  closely  follow  those 

Si 

for  the  core  enthalpy.  At  high  Mach  numbers  the  turbulent  results  are 

far  from  adiabatic,  and  the  theory  even  predicts  a  decrease  in  h  /H  if 

7  86 

transition  occurs  in  the  mixing  layer  at  Re^  less  than  10  (Fig.  20). 

Since  the  "fast  expansion"  distortion  of  the  initial  profiles  has  been 

neglected  in  these  calculations,  h  might  even  be  smaller  due  to  this 

(23)S 

effect,  if  the  laminar  results  carry  over  to  the  turbulent  case. 

Finally,  the  effect  of  Mach  and  Reynolds  numbers  on  the  base  pressure 
is  indicated  in  Figs.  21  and  22.  Generally  speaking,  the  range  of  P  /P 

D  00 
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is  about  .3  to  .6,  regardless  of  conditions.  One  notices  that  at  high 
Mach  numbers  the  laminar  and  turbulent  predictions  appear  to  agree 
within  about  15°/o.  If  the  theory  can  be  believed  under  these  condi¬ 
tions,  then  one  should  not  expect  transition  in  the  base  flow  of  a 
slender  re-entry  vehicle  to  be  accompanied  by  a  large  change  in  base 
pressure.  This  is  in  contrast  to  the  situation  for  Mach  3  adiabatic 
conditions  (shown  in  Fig.  22),  where  a  sudden  drop  in  base  pressure  is 
indicated  as  the  flow  changes  from  laminar  to  turbulent. 
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14.  CONCLUDING  REMARKS 


In  the  present  report  a  new  empirical  turbulent  mixing  and  base 
flow  model  has  been  constructed  which  hopefully  will  be  applicable 
slender  re-entry  vehicles.  A  plausible  method  has  been  found  to  extra¬ 
polate  turbulent  mixing  results  for  low  Mach  numbers  to  high  Mach  number 
flows  and  to  include  the  effects  of  non-similar  and  highly  cooled  turbu¬ 
lent  mixing  layer  development.  Reasonable  results  were  found  for  gross 
effects  such  as  the  base  pressure  and  rear  stagnation  enthalpy  for  a 
highly  cooled  10°half -angle  cone. 

The  asymptotic  mixing  layer  analysis  for  adiabatic  conditions  up  to 
about  Mach  3  can  probably  be  considered  reasonable,  since  the  velocity 
profile  data  appear  to  be  correlated.  Under  these  conditions  the  effect 
of  any  initial  boundary  layer  thickness  becomes  negligible,  since  the 
turbulent  profiles  develop  very  rapidly.  Application  of  the  present 
theory  to  any  situation  where  the  initial  profile  effects  may  become 
important  must  be  considered  conjecture  at  this  stage,  since  no  experi¬ 
mental  confirmation  is  presently  available.  This  includes  the  highly 
cooled  case  at  relatively  low  Mach  numbers  and  all  conditions  at  high 
Mach  numbers.  Hopefully,  the  present  treatment  of  non-similar  turbulent 
mixing  can  serve  as  a  guide  to  future  experiments  for  these  technically 
important  flow  cor, di  tions  . 
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APPENDIX  A 


Mavdew-Reed  Dividing  3treaa>line  Location 


Refers,  lag  to  Fig.  2,  it  may  be  seen  that  the  total  mass  flow  inside 
the  dividing  streamline  at  some  station  x  downstream  oi  the  nozzle  exit 
must  equal  the  mass  flow  through  the  nozzle.  Assuming  that  the  flow  is 
axisymmetr ic ,  the  mass  balance  yields 
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The  total  mass  flow  A  through  the  nozzle  was  calculated  by  assuming  an 
isentropic  expansion  from  the  reservoir  conditions  tabulated  by  Maydew 
and  Reed^  for  each  run. 


Since  the  experimental  porfile  data  u(r)  are  given  by  Maydew  and 

Reed  for  each  3tation  x,  the  integral  of  Eq.  (A- 1)  can  be  evaluated  as 

a  function  of  i *■  g  upper  limit  until  Eq.  (A-l)  is  satisfied.  This  o' 

course  requires  that  the  density  variation  be  known,  anf  this  was 

assumed  to  be  given  by  the  Crocco  integral  relation  for  a  perfect  gas, 

Eq.  (18).  The  dividing  streamline  velocity  is  then  given  by  u(r  „T). 

DSL 


An  estimate  of  the  accuracy  of  this  procedure  can  be  obtained  by 
assuming  that  the  flov  Inside  the  dividing  streamline  can  be  arbitrarily 
divided  into  a  '.rofile"  part  and  a  ''potential  core"  part.  If  most  of 
the  total  ft  ic  in  the  potential  core  (as  it  is  for  small  x) ,  then  a 
small  error  in  the  value  of  total  ft  will  have  a  large  effect  on  the 
value  of  For  example,  at  x  =  1.5"  a  5°/o  change  in  M  would  pro¬ 
duce  about  50  to  100°/o  change  in  u  for  the  experimental  conditions 

but  at  x  =  lj'1  a  5°/o  change  in  ft  would  produce  only  a  5  /o  to  1  ^  / o  in 
UDSL*  <*aCt  snovm  i°r  x  “  1*5  and  3.0  inches  therefore 

could  Le  considerably  in  error. 
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APPENDIX  E 


Eddy  Viscosity  Factor  k  of  Figure  3 


When  detailed  profile  data  are  available,  as  in  the  Maydew-Reed 
report,  the  factor  k  in  the  eddy  viscosity  expression  (Eq.  1)  can  be 
found  by  a  leas t -squares  fit  of  the  data  with  Eq,  ( 1 9) •  However,  much 
of  the  previous  literature  on  the  mixing  problem  does  rot  contain 
sufficiently  detailed  profile  data,  but  presents  only  the  final  result 
in  the  form  of  the  jet  spread  parameter  a,  A  sunmary  of  previous 
experimental  determinations  of  o  up  to  1962  is  given  by  Maydew  and  Reed, 
Since  a  depends  on  the  choice  of  profile  used  in  the  theory,  this  must 
also  be  specified. 

One  way  that  the  quantity  k  can  be  related  to  a  is  by  comparing 
the  derivatives  du*/d(y/x)  of  Eq.  (17)  and  (19)  at  some  selected  value 
of  u*.  For  example,  for  the  error  function  profile 


u*  =  1/2  fl  -  — 

\  f"“ 


(A- 2) 


we  match  the  slopes  at  u*  =  .5  to  get: 


(A-3) 


Thus  one  must  specify  the  profile  shape  function  g(cy/x)  and  the 
val ’e  of  u*  at  which  the  slopes  are  to  be  matched, 

.  learly  the  above  method  for  determining  k  may  not  be  very  accurate 
and  in  addition  depends  on  an  arbitrary  assumption  of  the  value  of  u* 
at  which  the  equr'-ion  is  evaluated,  This  procedure  therefore  will  not 
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be  used  In  this  report,  with  one  exception:  lhe  data  of  Zumwalt  at 
Mach  3  referred  to  by  Maydew  and  Reed  provide  an  additional  point  on 
Fig.  5  which  further  •'  mds  the  correlation. 

In  a  private  communication  with  Prof.  Zumwalt  at  the  University  of 
Oklahoma  it  was  found  for  the  error  function  profile  the  best  a  was 
loughly  23  to  30  at  Mach  2.9,  depending  on  how  the  initial  boundary 
layer  thickness  was  taken  into  account.  Eq.  (A-3)  was  evaluated  at 
u*  -  .5  using  this  information  to  give  the  data  points  of  Fig.  5 
attributed  to  Zumwalt. 


A  second  point  to  be  made  concerns  the  determination  of  k  by 
correlating  mixing  layer  profile  data  in  the  non-similar  growth  region. 
This  will  surely  be  a  problem  in  high  Mach  number  flows  if  the  relation 
of  Fig.  5  is  approximately  correct.  The  mixing  layer  experiment  should 
produce  either  velocity  c>r  density  profiles  (or  both)  as  a  function  of 
physical  x  and  y.  The  basic  equations  of  this  report  show  that  the 
theoretical  relation  for  u*(x,y,k)  is  given  implicitely  by  the  relations: 


y  -  y 


DSL 
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F*  =  F*(S*,u*) 


S*  =  S*(x* ,  Cf) 


Cc  =  Cc  tH  /H  .  h  /H  ,  R  ) 
f  f  we'ee  e 
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p  u  F  k  r  Jdx 
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F  =  F  (C  „ ,  H  /  H  ,  h  /  H  ,  •  i  ) 

w  w  i  ’  w  e  e  e  w 


r.  =  p(P.b) 

h  =  H  -  u'72 
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Thus  k  appears  explicitely  In  the  relation  for  x *  as  a  function  of  x. 
Although  somewhat  cumbersome,  these  relations  can  be  programmed  for 
a  computer  so  that  a  least  squares  determinat ion  of  k  from  velocity 
profile  da. a,  density  profile  data,  or  but 5:  car  1“  performed. 
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DIVIDING  STREAMLINE  VELOCITY  RATIO,  u/u( 


FIG.  1  MIXING  LAYER  REGIONS 


FIG.  2  GEOMETRY  OF  THE  JET  MIXING  EXPERIMENT 
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FIG.  3 

EXPERIMENTAL  DIVIDING 
STREAMLINE  VELOCITY 
CALCULATED  FROM  MASS 
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EDDY  VISCOSITY  FACTOR,  k  NORMAL  DISTANCE,  cry'/x 


FIG.  6 

EQUIVALENT  SKIN  FRICTION  COEFFICIENT 

C£  FOR  CONES  AND  WEDGES,  H  /H  =  0 
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FIG.  7 

EQUIVALENT  SKIN  FRICTION 
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FIG.  9 

TURBULENT  INITIAL  NTXING  LAYER 
PROFILE 
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FIG.  8 

EFFECT  OF  WALL  ENTHALPY  RATIO 
ON  C  FOR  R  =  I °7 


FIG.  10  DEVELOPMENT  OF  THE  VELOCITY  GRADIENT  PROFILES  WITH  STREAMWISE 
DISTANCE 
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ENTHALPY  FUNCTION,  W  ENTHALPY  FUNCTION 


VELOCITY  RATIO,  u/u#  VELOCITY  RATIO,  u/u. 
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FIG.  11  DEVELOPMENT  OF  THE  ENTHALPY  FUNCTION  W  PROFILES  WITH  STREAMWISE 
DISTANCE 
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FIG.  12  VELOCITY  ON  THE  DIVIDING  STREAMLINE 


FIG.  13  ENTHALPY  FUNCTION  W  Cv  THE  DIVIDING  STREAMLINE 
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MOMENTUM  THICKNESS  INTEGRAL,<1>  STREAMWiSE  DISTANCE  PARAMETER, 


RECIRCULATING  CORE  ENTHALPY, 


ENERGY  BALANCE  CONTROL  VOLUME 
(BOUNDARY  IS  INVISCID  ELO*  STREAMLINE) 


FIG.  16 

BASE  FLOW  MODEL 
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FIG.  17 

DIVIDING  STREAMLINE 
VELOCITY  BEFORE 
RECOMPRESSION 


FREF  STREAM  MACH  NUMBER.  M, 


FIG.  18 

RECIRCULATING  REGION 
ENTHALPY 


FREE  STREAM  MACH  NUMBER,  M, 
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BASE  PRESSURE  RATIO,  P,  /F,  NECK  CENTERLINE  ENTHALPY  RATIO 


FIG.  19  EFFECT  OF  MACH  NUMBER  ON 

CENTERLINE  STATIC  ENTHALPY 
AFTER  RECOMPRESSION 


FIG.  21  EFFECT  OF  MACH  NUMBER  ON 
BASE  PRESSURE 


FREE  STREAM  REYNOLDS  NUMBER, Re, 


FIG.  20  EFFECT  OF  REYNOLDS  NUMBER  ON 
CENTERLINE  STATIC  ENTHALPY 
AFTER  RECOMPRESSION 


FIG.  22  EFFECT  OF  REYNOLDS  NUMBER  ON 
BASE  PRESSURE 
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